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ABSTRACT 

LORAN-C  reliability  diagrams  produced  at  the  Defense 
Mapping  Agency  Hydrographic/Topographic  Center  depict  two 
types  of  data:  (a)  the  maximum  usable  groundwave  signal 
limit,  which  aids  the  LORAfJ-C  user  in  determining  which 
transmitters  provide  coverage  in  his  area  of  operation, 
and  (b)  the  predicted  uncertainty  of  a  LORAN-C  hyperbolic 
fix. 

Signal  limits  are  computed  using  Bremmer's  field  pre¬ 
diction  formula  (Ref.  1)  and  an  algorithm  that  predicts 
the  range  for  a  signal  of  predetermined  signal-to-r.oise 
ratio  propagating  along  an  electrically  inhomogeneous 
transmission  path  (Ref.  2).  Fix  uncertainty- predictions 
are  based  on  a  formula  relating  fix  uncertainty  to  (a) 
crossing  angle  between  lines  of  position,  (b)  system 
standard  deviation,  and  (e)  the  divergence  of  hyperbolic 
lines  of  position. 

Actual  range  and  fix  uncertainty  may  differ  from 
values  shown  on  reliability  diagrams,  depending  on  such 
factors  as  weather,  the  occurrence  of  geomagnetic  dis¬ 
turbances,  and  the  user's  direction  of  travel. 

Reliability  diagrams  currently  produced  show  signal 
limits  and  -fix  uncertainties  for  LORAN-C  chains  at  a-  scale 
of  1:5,000,000;  a  new  generation  of  reliability  diagrams 
could  show  data  at  a  reduced  scale  (1:10,000,000)  for  each 
LORAN-C  triad  (one  master  and  two  slave  transmitters), 
making  more  chain  and  transmitter  selection  information 
available  to  the  user. 


j  INTRODUCTION 

| 

i.  The  LORAK-C  Navigation  System 

i  ■  ... 

LORAN-C  is  a  pulsed,  low  frequency  (100  kHz)  elec- 
i  tronio  navigation  system  operated  by  the  U.S.  Coast  Guard, 
j  LORAN-C  provides  fix  data  to  vessels  operating  in  the 
I  northern,  northwestern,  and  central  Pacific,  the 
-j  Mediterranean,  the  northern  Atlantic,  and  the  U.S.  Coastal 
;  Confluence  Zone. 

j  LORAN-C  is  operable  in  either  of  two  modes:  (a) 

!  range-range  mode,  in  which  a  fix  is  defined  by  the  inter- 
.  section  of  two  circular  lines  of  positions  (LOP's)  defined 
j  by  the  arrival  times  of  LORAN-C  pulses  from  two  synchro- 
:  nized  transmitters,  and  (b)  hyperbolic  mode,  in  which  each 
'  LOP  is  hyperbolic  and  is  defined  by  the  difference  in 
|  arrival  times  of  pulses  from  two  synchronized  transmit- 
!  ters. 

Synchronization  of  LORAN-C  transmissions  is  main- 
.tained  by  system  area  monitors  which  record,  at  fixed 
i  positions,  time  differences  (TD's)  in  pulse  arrivals  from 
gieaoh  master-slave  transmitter  pair  operating  in  ths  area. 
*  If  TD's  measured  at  a  monitor  drift  excessively  from  the 
j  norm,  then  ths  monitor  directs  the  sieve  transmitter  to 
j  make  a  change  in  pulse  timing  to  compensate  for  the  drift. 


The  Ideas  expressed  In  this  paper  represent  the  opinions 
j  of  the  author  and  do  not  necessarily  reflect  official 
1  policies  of  the  DMA. 


A  LORAN-C  chain  consists  of  one  master  transmitter 
(which  initiates  the  transmission  of  pulses  from  other 
transmitters  in  the  chain)  and  several  secondary  (or 
slave)  transmitters;  each  LORAN-C  chain  is  assigned  a 
unique  group  repetition  interval  (GRI).  The  GRI  is  the 
length,  commonly  expressed  in  tens  of  microseconds,  of  the 
coded  sequence  of  pulses  _that?  comprise  the  LORAN-C  signal 
format. 


General  Description  of  the  LORAN-C  Reliability  Diagram 


Shown  in  Fig.  1  is  the  LORAN-C  reliability  diagram  for 
the  Gulf  of  Mexico  (GRI  7980).  (A  list  of  the  currently 
available  LORAN-C  reliability  and  coverage  diagrams  is  in 
Appendix  A.) 

Two  types  of  data,  depicted  on  a  Lambert  conformal 
conic  projection  of  land-sea  interfaces  at  a  scale  of 
1:5,000,000,  are  shown  in  the  diagram:  (a)  the  maximum 
usable  groundwave  signal  limits  for  signal-to-noise  ratios 
(SNR's)  of  1:3  and  1:10  and  (b)  predicted  fix  uncertainty. 

The  signal  limit  contours  define  the  areas  in  which 
the  LORAN-C  groundwave  signal  from  each  transmitter  in  the 
chain  is  of  sufficient  field  strength  to  be  detected  by 
either  commercially  available  receivers  (generally  cap¬ 
able  of  extracting  signal  at  SNR’s  of  at  least  1:3)  or 
military  receivers  (generally  capable  of  extracting  signal 
at  SNR's  of  at  least  1:10). 

Predicted  fix  uncertainty  contours  define  the  areas  in 
which  the  accuracy  of  the  LORAN-C  system,  operating  in 
hyperbolic  mode,  is  limited  by  random  errors. such  as  those 
caused  by  transmitter  instability  or  variable  radir.wave 
propagation  conditions.  Contours  representors  fix  uncer¬ 
tainties  of  500,  750,  and.  1500  ft.  are  shown  in  reli¬ 
ability  diagrams. 

The  information  shown  in  reliability  diagrams  aids  the 
LORAN-C  user  both  in  (a)  planning  chain  and  transmitter 
selection  for  a  voyage  and  (b)  making  enreute  changes  in 
chain  and  transmitter  selection  that  become  necessary  wnen 
a  transmitter's  signal  becomes  unusable  due  to  conditions 
such  as  transmitter  failure  or  change*  in  weather. 


Figure  1.  LORAN-C  reliability  diagram  for  the 
Gulf  of  Mexico  (GRI  7930).  (Signal  limits  for 
the  master  transmitter  and  for  slaves  W,  X, 
and  I  are  not  shown.) 


PREDICTING  THE  MAXIMUM  USABLE  GROUNDWAVE  SIGNAL  LIMIT 
Attenuation  of  the  LORAN-C  Signal 

The  LORAN-C  signal  loses  energy  as  it  is  transmitted 
along  the  earth's  surface;  this  loss  results  from  signal 
front  spreading,  energy  scattering  by  irregular  terrain, 
energy  absorption  by  the  earth  and  its  atmosphere,  etc. 
Factors  such  as  transmission  path  characteristics,  trans¬ 
mitted  power,  and  distance  traveled  by  the  signal  affect 
the  amount  of  energy  lost  by  that  signal.  Transmission 
path  character 1st ies  include  (a)  physical  properties,  such 
as  curvature  of  the  earth's  surface,  and  (b)  electrical 
properties  (ground  conductivity  and  permittivity,  atmo¬ 
spheric  refractivity)  which  vary  as  functions  of  weather 
(Ref.  3),  vegetation,  and  terrain  ruggedneas. 


LORAN-C  receivers  are  designed  to  sanple  the  LORAN-C 
'pulse  25  microseconds  (usee)  following  its  leading  edge; 
this  standard  sampling  point  (SSP)  occurs,  ideally,  at 
0.506  df  the  pulse's  peak  amplitude  (Ref.  4).  See  Fig.  2. 
For  purposes  of  predicting  field  strength,  the  signal 
level  of  the  LORAN-C  pulse  is  taken  to  be  the  root-mean- 
aquare  (rras)  .amplitude  of  a  continuous  wave  whose  ampli¬ 
tude  is  that  of  the  pulse  at  its  SSP  (Fig.  2). 

At  some  distance  along  its  transmission  path,  the 
groundwave  component  of  the  LORAN-C  pulse  loses  so  much 
energy  that  it  becomes  indistinguishable  from  the  ambient 
atmospheric  noise.  This  occurs  generally  at  a  SNR  of  1:3 
for  commercial  receivers  or  1:10  for  military  receivers. 
The  SNR  is  the  ratio  of  the  LORAN-C  signal's  field 
strength  at  0.128  of  its  peak  power  to  the  rms  field 
strength  of  the  ambient  atmospheric  noise,  as  determined 
from  CCIR  report  322  (Refs.  5,  6). 


/.o  L  Lp 


LOftAM-C- 


p‘'*t  o  0 

Amplitude. 

(per  cent)  -  i0  (t 


-*-701  *  SO  S, 
(rms  o(  Mai 

r  amplitude  of 

a  CO«t»mAo*s 
_ wave  uAtoic 

"  tkat  of  tie 
pul^e.  at 
its  SSP) 


Table  1  . 

Conductivity 
Type  of  Surface  (mhos/, meter) 

Sea  water  3  to  5 

Rich,  damp  soils, 
not  heavily  -  , 

leached  10*4  to  3xl0*4  15  to  30 

Dry  sandy  soils, 
or  heavily  , 

leached  areas  10"J 
Very  thin  soil,  ,  , 

over  rock  10  to  10”J 

Fresh  water,  - 

average  lake  10*J  to  3x10" 

Glacial  Ice  4x10-5 


Permittivity 
(relative  to 
free  space) 


10  to  20 


10"4  to  10"3 


5  to  10 


*3  to  3xl0’3  80 


4x10-5 


10  to  20 


r.eld 


(v«Hi/m)  |0-‘ 


Pulse  fength  ( p  sec) 

Figure  2.  The  LORAN-C  pulse  is  sampled  25  usee 
after  it  begins.  The  signal  level  of  the  pulse  is 
defined  as  0.707  of  0.506  of  the  pulse's  peak  am¬ 
plitude  (0.128  of  its  peak  power). 


Procedure  for  Predicting  Signal  Limits  Shown  on 
Reliability  Diagrams 

(Step  n 

A  number  of  azimuths  radiating  outward  from  each 
transmitter  in  a  LORAN-C  chain  are  chosen.  Each  azimuth 
represents  an  electrically  inhomogeneous  transmission 
path,  and  is  treated  as  a  "mixed  path"  consisting  of  a 
number  of  electrically  homogeneous  terrain  segments.  The 
electrical  properties  of  3ome  common  terrain  types  are 
shown  in  Table  1. 


Siy,»l  L*e*. t  (  r\.m  ) 

Figure  3.  The  signal  limit  for  each  type  of  ter¬ 
rain  segment  is  determined  using  field  strength 
curves  based  on  Bremmer's  formula.  In  this  ex¬ 
ample,  the  limit  for  a  signal  of  field  strength 
6x10'”  volt/meter  is  420  nautical  miles. 

(Step  3) 

Once  the  signal  limit  for  each  type  of  terrain  segment 
is  predicted,  the  signal  limit  of  the  inhomogeneous  trans¬ 
mission  path  can  be  predicted  using  the  method  derived 
below  (Ref.  2). 

Shown  in  Fig.  4  is  a  mixed-path  approximation  of  an 
inhomogeneous  transmission  path.  The  segment  lenghts  are 
labeled  L(i,J);  the  segments  are  numbered,  outward  from 

the  transmitter,  i=  1 ,  i=2 . i=n  for  terrain  types 

J-1,  j-2,  ...  ,  Jam. 

Having  predicted J;he.signal  limits  S(j)  for  each  ter- 
rain^type  and  S  (3ea)  for  an  all  sea  water  path,  a  sea- 
equivalency  Tactor  TTJ)  is  defined: 


F(J)=S(sea)/S(J) 


(Eq.  1) 


Each  terrain  segment  of  length  L(i,J)  can  be  converted  to 
a  sea-equivalent  segment  of  length  G(i)  using  the  relation 


G(l)sF(J)  X  L(i, j) 


(Eq.  2 


(The  length  G(i)  is  the  distance  a  signal  would  travel 
over  sea  water  before  attenuating  to  the  level  it  would 
attenuate  to  by  traveling  a  distance  L(i,j)  over  terrain 
of  type  J.) 
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(Taken  from  Ref.  7) 

Tna  division  of  an  azimuth  into  homogeneous  terrain 
segments  is  based  on  maps  of  either  conductivity 
boundaries  or  land-sea  Interfaces.  (Shown  in  Ref.  8  is  a 
map  of  conductivities  for  the  United  States.) 

(Step  2) 

Ones  the  average  ambient  atmospheric  noise  level  for 
the  area  surrounding  a  transmitter  is  established,  the 
field  strengths  of  the  LORAN-C  signal  at  SNR's  of  1:3  and 
1:10  are  calculated.  Ualng  field  strength  curves,  suoh  as 
the  one  in  Fig.  3,  the  signal  limit  of  eaoh  type  of  terrain 
segment  is  detsrminsd. 

Field  strength  curves  are  computed  using  Breaker's  far 
field  prediction  formule  (Appendix  B). 


Figure  4.  Mixed-path  approximation  of  an  inhomo¬ 
geneous  transmission  path. 


Applying  Eq.  2  eonverts  a  mixed-terrain  transmission 
path  to  an  all  sea  water  transmission  path.  The  sea  water 
signal  limit  S(sea)  is  known,  so  the  mixed-path  signal 
limit  can  bs  calculated  by  converting  all  sea-equivalent 
segments  up  to  S(ssa)  back  to  terrain  segments,  then  sum- 
sing  their  lengths: 


2  C(i)  *  S(ses) 
1*1 


(Eq.  3) 


where  the  A'th  ses-equivslent  segment  extends  to  S(sea). 

If  the  A'th  segment  "straddles"  S(ses),  in  other 
words, 


52  0(1>  >  «*•■> 


(Eq-  «>  ! 


and 

« 


A*  I 

XC(i) 

i«l 


S(sea) 


(Eq.  5) 


then  Eq.  3  must  be  modified. 

ktJh*5/Par>  of  th*  A'th  3«*-«quivalent  segment  that  lies 
within  S(sca)  is  the  residual 


A-1 

RES  s  S(sea)  -  ^  G(i) 
i>1 


The  aixsd-path  signal  limit  S(aixed)  is  predicted  by 
converting  the  sea-equivalent  segments  Gill  through  GCA-1 ) 
and  RES  back  to  their  original  lengths  Mi,  j)  and  r.  and 
|  summing: 


where 


S( mixed)  = 


riRES/F< j) 


A-1 

♦  r 

ill 


(P(J)A  i  F(j)  for  iiA). 

Thus, 

S(m.*<ul)'-  Z  L(i,iV  F(j)lil,}) 

«j>a 

i  As  an  example,  refer  to  Tables  2"  and  3.  In  Table  2  are 
predicted  signal  limits  and  sea-equivalency  factors  for 
three  terrain  types;  in  Table  3  are  sea-equivalent  lengths 
C(i)  calculated  for  each  of  five  segments  in  a  mixed- 
terrain  transmission  path.  (Note  that  the  fourth  terrain 
segment  straddles  S(sea),  hence  A=4.)  Using  Eq.  10,  the 
mixed-path  signal  limit  is  calculated  to  be  353  nautical 
miles  (n.-i.).  This  technique  for  predicting  mixed-path 
signal  limits  is  represented  graphically  in  Fig.  5. 
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Table  2 

Sea-Equi¬ 
valency 
factor  F(,1) 

1.1 

4.5 

1.0 


Signal 
Limit  S(J) 
(n.ml. ) 

810 

200 

900 


(Eq.  6). 
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Figure  5.  Prediction  of  mixed-path  signal  limits 
shown  on  reliability  diagrams:  (1)  approximate 
inhomogeneous  transmission  path  as  a  mixed-path, 

(2)  convert  terrain  segments  to  sea-equivalent 
segments  and  determine  length  of  residual  RES,  and 

(3)  convert  sea-equivalent  segments  and  RES  back 
to  terrain  segments  (Ref.  2). 

(Step  4) 

For  each  azimuth  (transmission  path)  radiating  outward 
from  a  given  transmitter,  predicted  3ignal  liait3  for 
SNR's  of  1:3  and  1:10  are  plotted,  then  contoured. 


j  Other  Methods  of  Predicting  Mixed-Path  Signal  Limits 

t 

Using  Eq.  10  to  predict  signal  limits  requires  much 
less  computer  tine  than  other  known  methods,  such  as 
Millington'3  (Refs.  9-12).  Millington's  method  of  pre¬ 
dicting  the  field  strength  of  a  groundwave  transmitted 
over  a  mixed-terrain  path  is  given  by 


T  E,(d,+dJ 
Erdj+d2+...+dJ  ‘1  E,  (d,  )_UL_1  ■ 


^  V-dn> 


lv  1' 


w 


w,'l 

E  , (d  ,+d  ) 
n-1'  n-1  n' 

E  ,(d  )  £,(d  +d  . .+d0) ! 

n-1'  ir  1'  n  n-1  2'J 


En(d1+d2*>..+dn) 


En(dl+--'+dn-l-) 


Ei(dj+d2+ .. +dn) 


(Eq.  11) 


where  E,(d.)  is  the  field  strength  of  a  groundwave  trans-  j 
mitted  a  distance  d  over  a  terrain  of  type  j. 

To  predict  a  signal  limit  using  Millington's  method, 
Eq.  11  must  be  solved  for  a  range  of  transmitter-to-  | 
receiver  distances  to  determine  at  what  distance  the  sig-  . 
nal  attenuates  to  its  1:3  or  1:10  SNR  level.  j 

Signal  limits  predicted  using  Eq.  10  were  within  14% 
of  the  limits  predicted  using  Millington's  method  in  17 
of  the  20  representative  mixed-paths  tested.  Eq.  10  pro¬ 
duced  signal  limit  predictions  that  were,  on  the  average, 
11%  more  pessimistic  than  those  produced  by  Millington's 
method. .  j 


Note:  -  S(sea) 

*  3(3). 

Tabla  3 

Terrain 
Segment  1 

Terrain 
Type  J 

Sea-Equivalent 
Segment  Length  Length  G(l) 

L ( 1 . 1 )  (n.ml.)  (n.ai.) 

1 

2 

50 

22  5 

2 

3 

100 

100 

3 

T 

100 

110 

4 

2 

500 

2250 

5 

3 

600 

600 

Eq.  10  does  not  take  into  account  the  fact  that  a 
signal's  attenuation  rate  (in  decibels/ 1 00 -n. mis-,  -for- in- 
■stance)  -  decreases  with  distance  from  the  transmitter  at 
different  rates  for  different  kinds  of  terrain.--  The  ef- 
feot  this  has  on  computing  sea-equivalent  segments  is 
shown  in  Table  4:  a  terrain  segment  has  a  sea-equivalent 
|  segment  whose  length  varies  with  distance  from  the  trans- 
!  aitter.  One  method  of  predicting  attenuation  over  a  mixed 
path  that  accounts  for  the  effects  of  the  variable  at¬ 
tenuation  rate  is  described  below  (Ref.  14). 

Consider  the  terrain  segments  described  in  Tables  2 
and  3-  Using  the  set  of  field  strength  curves  m  Fig.  6, 

the  attenuation  (in  decibels, -dB)  of  a  signal  propagating 
over  terrain  type  2  between  0  and  50  n.ml.  from  lta  trans¬ 
mitter  is  determined.  Next,  the  attenuation  of  a  signal 
propagating  between  50  and  150  n.ml.  Is  determined,  and  so 
on.  Tha  total  attanuation  over  a  distance  d  ia  tha  sum  of 
the  attenuation  valuaa  (in  dB)  of  all  segments  up  to  dis¬ 
tance  d. 


Note:  RES  «  900-(225*100-»I10)  *  465  n.ml 


Table  4 


Attenuation  of 
Signal  Transmitted 
Across  Terrain  Sea-Squivalent 

Terrain  Segment  Segment _ •  Segment _ ” 


1  to  101  1cm  41  dS 

400  to  500  3.5 

1000  to  1100  2.5 


1  to  115  km 
400  to  520 
1000  to  1150 


Figure  7.  Fix  undertainty  due  to  LOP  instability 
over  distances  a  and  b  is  exaggerated  by  the  small 
crossing  angle  A. 


*  assuming  permittivity  (relative  to  free  space) 
is  15  and  conductivity  is  10-'3  mho/meter. 

••  sea-equivalent  segment  length  determined  for 
the  attenuation  given  in  column  two. 

Data  taken  from  Ref.  13. 


"D  utanct  frt *«  Trar<i*niii  er 

( n-mi  ) 


Figure  6..  The  total  attenuation  over  the  mixed 
path  from  0  to  d  is  A*B+C*D  decibels. 


PREDICTION  OF  LCHAN-C  FIX  UNCERTAINTY 
Definition  of  Fix  Uncertainty 

A  LOR AN -C  position  fix  is  defined  by  the  intersection 
of  two'  hyperbolic  LOP's.  Due  to  factors  such  as  trans¬ 
mitter  instability  and  variable  propagation  conditions, 
i  the  geographic  location  of  each  LOP  changes  with  time; 
this  variation  is  essentially  random  and  is  assumed  to 
occur  as  a  normal  distribution.  The  location  of  the  in¬ 
tersection  of  two  LOP's  is  thus  assumed  to  vary  as  a  two- 
dimensional  normal  distribution. 

P'x  uncertainty  contours  depicted  on  L0RAN-C  reli- 
eoillcy  diagrams  show  area  in  which  a  fix  has  a  951  chance 
of  being  within  1500,  750,  or  500  ft.  of  the  center  of  this 
two-dimensional  normal  distribution,  barring  the  occur¬ 
rence  of  non-random  errors.  (The  center  of  the  distribu¬ 
tion  may  or  may  not  be  associated  with  a  distinct  geodetic 
position. ) 

As  an  example,  consider  a  LGRAN-C  receiver,  located  at 
a  stationary  position  within  the  area  defined  by  the  500 
ft.  uncertainty  contour,  monitoring  TD's  for  several 
|  weeks.  95$  of  the  position  fixe3  defined  by  these  TD's 
|  will  be  within  500  ft.  of  the  stationary  receiver. 

I  Fix  uncertainty  due  to  the  random  instability  in  LOP 
I  location  Increases  as  the  crossing  angle  between  LOP's 
|  decreases  (Fig.  7).  In  addition,  fix  uncertainty  in¬ 
creases  as  hypcrbolio  LOP's -diverge.  For  example,  a  TD 
I  error  of  0.1  usee  along  the  baseline  between  a  master  and  a 
,  slave  transmitter  corresponds  to  a  distanoe  error  of  about 
!  50  ft.;  near  the  baseline  extension,  where  LOP's  are 
|  highly  divergent,  the  same  TD  error  represents  a  distance 
,  error  of  several  hundred  feet  (Fig.  8). 

1 


Figure  8.  Hyperbolic  LOP's  (labeled  in  micro¬ 
seconds  of  time  difference)  diverging  away  from 
the  baseline. 


Thus,  fix  uncertainty  contours  show  areas  in  which  the 
i  accuracy  of  the  L0RAN-C  system  is  limited  by  the  random 
|  instability  of  hyperbolic  position  lines.  Fix  uncertainty 
due  to  this  instability  is  amplified  (a)  where  crossing 
i  angles  between  LOP's  are  small  and  (b)  wnere  LOP’ 

I  diverge. 

Fix  uncertainty  must  not  be  confused  with  "absolute 
accuracy"  -  the  accuracy  with  which  a  LORAN-C  fix  can  be 
associated  with  a  distinct  geodetic  position.  The 
absolute  accuracy  of  a  LORAN-C  fix  is  a  function  of  both 
(a)  fix  uncertainty  and  (b)  the  ability  te  correctly  pre¬ 
dict  the  time  required  for  a  LORAN-C  signal  to  propagate 
from  its  transmitter  to  any  given  geodetic  position. 

Fix  uncertainty  is  similar  to  repeatable  accuracy  (the 
accuracy  with  which  a  vessel  can  return  to  a  previous 
position  using  LORAN-C).  In  an  area  of  500  ft.  repeatable 
accuracy,  a  vessel,  using  the  same  LORAN-C  TD  coordinates, 

<  can  return  to  the  same  position  repeatedly  to  within  500 
ft.,  most  of  the  time. 

In  an  area  of  500  ft.  fix  uncertainty,  a  vessel  will  te 
j  able  to  position  itself  to  within  500  ft.  of  a  buoy  using 
•  the  TD  coordinates  of  that  buoy,  provided  these  TD  coor¬ 
dinates  have  been  established  by  monitoring  TD's  at  the 
buoy  for  a  sufficient  length  of  time. 


i  Derivation  of  an  Approximate  Method  for  Predicting  Fix 
j  Uncertainty 

Fix  uncertainty  data  depicted  on  LORAN-C  reliability 
diagrams  are  computed  using  a  formula  derived  by  Trow  and 
Jessell  (Ref.  15).  A  similar  formula  is  derived  by 
Sltterly  (Ref.  16).  This  formula  defines  the  "95  percent 
radial  error"  -  the  radius  of  the  circle  containing  about 
95$  of  all  fixes  associated  with  a  given  TD  pair. 

Shown  in  Fig.  9  are  two  L0Prs  intersecting  at  an  angle 
A  and  displaced  from  their  average  positions  by  random 
errors  e(1)  and  e(2).  The  major  diagonal  d  of  the  paral¬ 
lelogram  thus  formed  is  the  position  fix  error  associated 
with  random  errors  e  (1)  and  e(2).  d  can  be  expressed, 
according  to  the  Law  of  Cosines,  as 
2  2  2 

d  «  P  ♦  q‘  ♦  2pq  cos A  (E*.  12) 


(p  and  q  are  tha  lengths  of  tha  parallalograaa  sides) 


By  the  definition  of  aine 


c(A 

o  -  t(i) 

(Eq.  13). 

P  *  Sin  A 

'  V  sin  A 

Combining  Eqs. 

12  and  13  gives 

dl-  L&Q  J 

)*,.  zet'itd)  cos  A 

(Eq.  14) 

6  Vs. n  A/ 

*  Sin  A  1  S'OlA 

Eq.  14  relates  position  fix  error  d  to 
variations  that  occur  in  LOP  location  and  (b) 

(a)  random 
LOP  crossing 

r(2).  Because  or  instability,  the  LOP  at  R  may  be  dis¬ 
placed  a  distance  DISP  to  R';  DISP  can  be  approximately 
expressed  in  terms  of  a  variation  dr  in  ray  path  lengths 
r(l)  and  r(2): 


DISP 


dr  /  B 


angle.  The  magnitudes  of  roughly  95<  of  these  position 
fix  errors  are  smaller  than  twice  their  root  mean  square 
(twice  their  standard  deviation)  2drms.  (This  approxima¬ 
tion  is  made  assuming  that  the  fix  error  d  is  normally 
distributed,  which  it  is  not;  the  probability  that  d  is 
smaller  than  2drms  varies  with  the  crossing  angle.)  Thus, 
fix  uncertainty  is  approximated  as 

r«d,.l  srrof  X  2Jr»(  =  zViTZdT  (Eq.  15). 


(Eq.  20) 

The  hyperbolic  position  line  at  R  is  defined  by 

r(l)  -  r(2)  s  constant  (Eq.  21). 

So  the  displaced  LOP- at  R*  is  defined  by 

(r( 1 )*dr)-(r(2)-dr)  *  r(1 )-r(2)+2dr  =  constant  (Eq.  22) 

The  distance  2dr  can  be  expressed  as  a  time  difference 
error  dt: 


dt  =  2dr/c 


(Eq.  23) 


where  e  is  the  speed  of  light  through  the  atmosphere. 
Combining  Eqs.  20  and  23  gives 
'/%  dt  •  C 

D,SP=  “slnB  (Eq*2U)- 

Letting  DISP  be  a  random  error  in  LOP  location  (such  as 
e(1)  or  e(2)  in  Fig.  9)  the  standard  deviation  sd  becomes 

sat- 1 


s>n  S 


(Eq.  25) 


Figure  9.  d  is  the  position  fix  error  associated 
with  random  errors  e(1)  and  e(2). 


Substituting  from  Eq.  14,  Eq.  15  becomes 


where  sdt  is  the  standard  deviation  of  random  time  dif¬ 
ference  errors  (also  called  the  "system  standard  devia¬ 
tion")-. 

Thus,  assuming  no  correlation  in  random  errors  (C=o; . 
Eq.  19  (which  approximates  fix  uncertainty)  reduces  to 


Zivmt  « 


Silt  -  c 


A  jisin  86))  +  (tin  (Eq-  26). 

(1/sinB(1)  and  1/sinB(2)  are  often  called  "expansion  fac¬ 
tor S"  or  "position  line  divergence  factors".) 


(Eq.  16). 


I 


Recalling  that  the  standard  deviation  of  a  normal  distri¬ 
bution  of  errors  is 


(Eq.  17) 


and_ letting  the  "correlation  coefficient"  be  defined  as 


C» 


gefrjeffl 


(Eq.  18), 


Eq.  16  becomes 


Figure  10.  The  displacement  DISP  of  an  LOP  can  be 
expressed  approximately  in  terms  of  the  deviations 
dr  in  raypath  lengths  r(1)  and  rj2).  (Taken  from 
Ref.  15.) 

Fix  Uncertainty  Contours 


Zirms  *  Shbr(s«K''S1*  sift1*  lsAttsAXi)Ce.*Ajl  (Eq.  19). 

The  correlation  coefficient  C  varies  between  -1  and  1 , 
indicating  the  degree  tor  which  random  errors  e(D  and  e(2) 
ere  related  to  one  another  by  some  common  causal  mecha- 
niam. 

The  standard  deviations  add)  and  sd(2)  in  LOP  loca¬ 
tion  should  be  expressed  in  a  way  that  reflects  the 
effects  of  LOP  divergence.  Shown  in  Fig.  10  srs  two  ray 
path  distances  r(1)  end  r(2)  between  eseh  of  two  transmit¬ 
ters  and  a  LORAN-C  reoelvtr  (located  at  R).  Also  shown  is 
the  line  segment  tangent  to  the  hyperbollo  LOP  at  R.  This 
assent  bisects  the  angle  2B  between  ray  paths  r(1)  and 


|  .For  a  LORAN-C  triad,  the  crossing  angle  A  is  the  sum 
(or  difference)  of  angles  B(1.)and  B(2),  as  shown  in  Fig. 

I  11.  Thus,  for  a  given  location  relative  to  the  LORAN-C 
I  triad,  the  fix  uncertainty  2dros  can  be  computed  using  Eq. 
'26. 

Fix  uncertainty  contours  shown  on  LORAN-C  reliability 
1  diagrams  are  derived  by  first  solving  Eq.  26  for  each  node 
:  in  a  spheroidal  grid  of  latitudes  and  longitudes  in  the 
vicinity  of  the  LORAN-C  triad.  Next,  an  interpolation 
1  scheme  is  used  to  determine  the  geographic  coordinates  for 
contours  of  1500,  750,  and  500  ft.  fix  uncertainty.  These 
coordinates  are  plotted,  using  a  Lambert  conformal  conic 
projeotion,  for  eaoh  triad  in  a  LORAN-C  chain;  the 
resulting  fix  uncertainty  plots  srs  eomposltsd  for  the 
entire  LORAN-C  ohsln. 

It  is  assumed  in  these  computations  that  the  standard 
|  deviation  of  TD  arrors  sdt  is  0.1  useo. 
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Figure  11.  Crossing  angle  A  is  the  sum  (or  dif¬ 
ference)  in  bisecting  angles  B(1)  and  B(2). 
(Taken  from  Ref.  15.) 
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A  More  Exact  Method  of  Predicting  Fix  Uncertainty 


The  error  d  due  to  random  variations  in  LOP  location 
does  not  occur  as  a  normal  distribution.  The  actual  dis¬ 
tribution  of  errors  is  described  by  the  relation 


qsJb  radial  error  *  kdi-ms 


:JTZd' 


where  K  varies  between  about  1.73  and  1.96  as  a  function  of 
the  crossing  angle  and  the  ratio  3d(1)/sd(2).  Approxi¬ 
mating  the  95J  radial  error  as  2.00  dras  thus  produces  a 
pessimistic  prediction  of  fix  uncertainty. 

In  an  approach  desribed  by  Sltterly  (Ref.  16)  and  in¬ 
vestigated  in  detail  by  Hiraiwa  (Refs.  17,  18)  the  coef¬ 
ficient  K  is  computed;  this  two-part  process  is  described 
below. 

(Step  1) 

Shown  in  Fig.  12  are  two  LOP's  intersecting  at  an 
angle  A.  The  probability  thar a  position  fix, 'associated 
with  the  intersection  of  L0P(1)  and  L0P(2),  will  fall 
within  the  very  thin  parallelogram  PQ  is  given  by 
aft  v(i1 

P'.k.u-Ay  =  — Jwsit j«p(f&)4v  (Eq' 28K 

uM  v(>) 

sd(1)  and  sd(2)  are  the  standard  deviations  of  L0P(1) 
and  L0P(2)  from  their  average  locations,  du  is  the  thick¬ 
ness  of  PQ  in  the  direction  perpendicular  to  L0P(1). 
Lengths  v(1)  and  v(2)  are  perpendicular  to  L0P(2)  and  de¬ 
fine  the  points  P  and  Q. '  As  shown  in  Fig.  12,  ”v(1)  and 
v(2J  .are  defined,  in  terms  of  perpendicular  distance  u 
from  L0P( 1 ) ,  as 


u  cosA  *  VZ1-  u*  sin  A 


(Eq.  29) 


The  probability  that  a  position  fix  lies  within  a 
circle  of  radius  Z  can  thus  be  computed  using  Eqs.  28-30, 
numerically  Integrating  over  the  Interval  u(1)*0,  u(2)*Z, 
and  multiplying  the  results  by  2.  This  is  done  for  various 
values  of  Z,  sd(1)/sd(2),  and  crossing  angle  A  to  deter¬ 
mine  which  combinations  of  these  values  produce  a  proba¬ 
bility  of  951. 


Figure  12.  The  probability  that  a  position  fix, 
defined  by  L0P(1)  and  L0P(2),  will  fall  within  the 
very  thin  parallelogram  PQ  is  given  by  Eq.  28 
where  v(1)  and  v(2)  are  defined  as  above. 
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Figure  13.  Graph  of  K  as  a  function  of  crossing 
angle  and  the  ratio  sd(1)/sd(2).  (Taken  from  Ref. 

17.) 

Having  determined  K,  the  "true"  951  radial  error  is 
computed  as 


Kdri»s 


'4  k  sJt-c 
Sir.  A 
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(See  Eq.  26). 

A  comparison  of  fix  uncertainty  values  computed  using 
Kdrms  with  those  computed  using  2drms  is  shown  in  Fig.  14. 

Reliability  diagrams  produced  in  the  future  by  DMA  may 
depict  fix  uncertainty  data  computed  using  Kdrms. 


(Step  2) 

Using  the  same  steps  by  which  Eq.  19  is  derived,  Eq.  27 
becomes  (for  C*0) 

zt  95?*  rsWi'el  ernjr  a  “\J sA(>)*+  *A(i)1  (Eq.  31). 

>q.  31  Is  solved  for  K  using  ths  various  values  of  Z, 
sd(1)/sd(2),  and  A  determined  in  Step  1.  Shown  in  Fig.  13 
is  s  graph  of  X  as  a  function  of  ad(1)/sd(2)  and  A. 


PRESENTATION  OF  DATA 
Reliability  Diagrams 

Fix  unesrtalnty  and  signal  limit  data  are  plotted  for 
each  L0RAN-C  chain  on  a  1:5,000,000  Lambert  conformal 
eonie  projection  of  land-sea  interfaces.  No  signal  limits 
are  shown  outside  the  1500  ft.  fix  uncertainty  contour. 


Figure  14.  Comparison  of  fix  uncertainty  values 
computed  as  Kdrms  with  those  computed  as  2drms; 
note  that  the  difference  is  most  significant  for 
small  crossing  angles. 


i  DIFFICULTIES  AND  IMPOSSIBILITIES  IN  PREDICTING  SIGNAL 
LIMITS  AND  FIX  UNCERTAINTIES 

i  Actual  signal  limits  and  fix  uncertainties  often  dlf- 
j  fer  from  predicted  values,  due  to  a  variety  of  causes, 

;  including: 

(a)  our  inability  to  predict  either  the  occurrence,  or 
the  effects  on  signal  propagation,  of  non-periodic 
phenomena  such  as  thunderstorms  and  geomagnetic 
disturbances; 

I  (b)  our  inability  to  accurately  predict  the  effects  of 

j  .  periodic  phenomena,  such  as  seasonal  climate 

I  changes,  on  signal  propagation; 

1  (c)  Inaccuracies  in  predictions  of  noise  level, 

|  groundwave  field  strength,  and  fix  uncertainty  due 

i  to  assumptions  made  in  the  preditlon  models  (for 

1  example,  the  system  standard  deviation  will  differ 

from  the  assumed  0.1  usee,  depending  on  the  pro¬ 
pagation  conditions  that  exist  at  a  given  moment); 
!  (d)  limitations  in  availability  and  accuracy  of  data, 

such  as  ground  conductivities,  used  in  the  predic¬ 
tion  models. 


i  Coverage  Diagrams  i 

I  1 

j  LORAN-C  coverage  diagrams  show  contours  of  skywave  and 
I  groundwave  coverage.  Groundwave  coverage  contours  define 
j  areas  in  which  these  conditions  exist: 

j  (a)  predicted  fix  uncertainty  is  1500  ft.  or  less,  and 
I  (b)  predicted  groundwave  3ignal  strength  exceeds  a  1:3 

SNR.  I 

Triad  Diagrams  I 

It  is  often  impossible,  using  currently  available  '■ 
reliability  diagrams,  to  determine  fix  uncertainty  for  a 
given  triad  within  a  LORAN-C  chain.  For  example,  the  dark  ; 
solid  line  in  Fig.  15  is  the  750  ft.  fix  uncertainty  con¬ 
tour  for  the  chain  HXYZ.  Using  this  chain  diagram,  a  user  ; 
operating  in  the  shaded  area,  using  TD's  from  the  MX!  j 
triad,  cannot  determine  fix  uncertainty.  a  triad  diagram  i 
for  transmitters  M,  X,  and  Y  reveals  the  location  of  the  , 

750  ft.  contour  (dashed  line  in  Fig.  15).  j 

Four  triad  diagrams,  each  at  a  scale  of  1:10,000,000,  j 
can  be  printed  on  the  same  size  sheet  as  one  of  the  cur-  : 
rently  available  reliability  diagrams.  j 


Figure  15.  750  ft.  fix  uncertainty  contour  (dark 
solid  line)  for  the  MXYZ  chain.  A  user  operating 
in  the  shaded  area,  using  70's  from  the  MX 7  triad, 
cannot  determine  fix  uncertainty  using  this  dia¬ 
gram. 


Absolute  Accuracy  J 

In  areas  whera  LORAN-C  data  are  controlled  by  high 
quality  data  from  systems  such  as  NNSS  (Transit  satellite) 
or  Autotape,  accuracy  information  could  be  depicted  in 
reliability  diagrams,  possibly  as  shading  in  various  tones 
of  gray.  I 


In  addition,  variables  such  as  the  user's  direction  of 
travel  affect  signal  limits  and  fix  uncertainties  observed 
by  the  user  (a  LORAN-C  receiver  may  "lock  on"  to  a  signal 
close  to  the  transmitter  and  track  that  signal  far  beyond 
its  predicted  limit). 

In  an  effort  to  make  reliability  diagrams  as  accurate 
and  useful  as  possible,  I  make  the  following  recommenda¬ 
tions: 

(a)  use  the  method  described  in  Fig.  6  to  predict 
groundwave  signal  limits; 

(b)  use  Hiraiwa's  method  (Eq.  32)  to  predict  fix  u: 
certainty; 

(c)  present  fix  uncertainty  and  signal  limit  data  for 
single  triads,  rather  than  for  entire  chains; 

(d)  where  data  are  available,  show  absolute  accuracy 
information  on  reliability  diagrams. 


APPENDIX  A.  CURRENTLY  AVAILABLE  LORAN-C  COVERAGE  AND 
RELIABILITY  DIAGRAMS 

DMA  Stock 

Miimhan  Area  (GRI) 


WOBZP5130 

LORAN-C  Coverage  Diagram 

WOBZP5131 

LORAN-A  Coverage  Diagram 

Reliability  Diagrams: 

2L0RX5592 

Mediterranean  Sea  (7990) 

ZLORX5593 

Norwegian  Sea  (7970) 

ZLORX5595 

North  Pacific  (9990) 

ZLORX5596 

Central  Pacific  (1990) 

ZL0RX5597 

Northwest  Pacific  (9970) 

ZLORX5598 

North  Atlantic  (7930) 

ZL0RX5600 

Gulf  of  Alaska  (7960) 

ZLORX56 01 

Canadian  West  Coast  (5990) 

ZL0RX5602 

West  Coast,  U.S.A.  (9910) 

ZL0RX5603 

Southeast  U.S.A.  (7980) 

ZL0RX5601 

Northeast  U.S.A.  99960) 

ZL0RX5605 

Great  Lakes  (8970) 

!  These  diagrams  are  available  through  agents  of  the  DMA 
|  Office  of  Distribution  Services;  addresses  of  these  agents 
are  listed  in  the  Defense  Manning  Aveney  Catalogue  of 
Haps,  Charts,  and  Related  Products!  Publication  1-ii-L,  DMA 
stock  number  CATD1NL). 


APPENDIX  B.  BREMMER'S  FAR  FIELD  PREDICTION  FORMULA 


k  i  =  permittivity  of  the  earth*  The  relative  permit-  ' 
tivlty  C,/f.a  is  about  15  for  land,  and  80  for  sea.  ; 

6  -  conductivity  of  the  earth  (mho/neter).  and  £ 

vary  a3  functions  of  weather,  soil  conditions,  etc. 
Conductivity  is  a  measure  of  the  electrical  "absorbancy"  ; 
of  a  medium. 


A»  =  permeability  of  free  space 
henry /meter. 


1.26  x  10" 


Brcmmer’s  formula  (Refs.  1,  9,  10,  19,  20,  21,  22)  is  • 
used  to  predict  the  vertical  field  of  a  groundwave  prop-1 
agating  along  a  smooth,  spherical,  electrically  homo-' 
geneous  earth.  The  strength  of  this  field  (in 
volts/meter)  is  calculated  as 

)Er|=  (E(R)2  +  E(I)2)*  (Eq.  33) 


where  E(R)  and  E(I)  are  the  real  and  imaginary  components 
c.'  tho  ground*:.. ve*  fiolil. 

Kcj*  3 4  ^r  ru:'r^ula ,  below  which  are  Jcfirwcl 

■the  terras  used  in  the  formula. 


tr- 


‘frre.oJd 

^  e’((k. TZl  *  ¥ 

n'o  -  W 


(Eq.  34) 


Er  a  vertical  electrical  field  (volts/meter),  both  trans¬ 
mitter  and  receiver  on  the  ground. 

i.  VT 

lol  =  dipole  moment  of  source  (amp-meters) .  Computations 
can  be  made  assuming  lol  =  1,  then  the  results  scaled  Ujl 
according  to  the  relationship  (Ref.  23)  Er  =  0.3VP 
volts/meter  at  1  ks  (0.54  n.mi.)  from  the  transmitter. 
P  is  0.128  of  the  peak  transmitted  power  in  kilowatts. 

Kj  «  ‘tf.'/ii  -  wave  number  of  the  atmosphere 
1  (radians/meter ). 

w  *  2fr  f  =  radial  frequency  of  transmitted  signal. 

f  *  linear  frequency  of  the  transmitted  signal.  No  core 
than  1?  of  the  radiated  energy  in  a  LC5.AN-C  signal  lies 
outside  the  90-110  kHz  band  (Ref.  2d). 

c.  =  speed  of  light  in  free  space  s  2.99792458  x  108 
meters/second. 

*/S  s  index  of  refraction  of  the  atmosphere  =  1.000338 
at  the  earth's  surface  for  a  standard  atmosphere.  Y*ET 
varies  as  a  function  of  weather. 

r  permittivity  of  the  atmosphere. 

-12 

€»  *  permittivity  of  free  space  *  8.85  x  10 

farad/raeter.  Permittivity  is  the  measure  of  electrical 
"transparency"  of  a  propagation  medium. 

d  =  distance  travelled  by  signal.  Eq.  34  is  useful  for  d 
>80  km  for  L0RAN-C  frequency. 

*  effective  radius  of  earth  =  4/3  of  radius  of  a 
spherical  earth.  ae  is  used  in  Eq.  34  to  compensate,  in 
part,  for  the  effects  of  an  inhomogeneous  atmosphere. 

each  X,  term  has  a  real  and  an  imaginary  component  and  is 
approximated  as  the  root  of  Ricattl's  differential 
equation  Jt/Jr,  -ll'ti'l*  0  •  A  detailed  explanation  of  how 
Xt  Is  determined  is  given  in  Ref.  22  . 

kg  b  ^6.,+  .* L  s  wave  number  of  the  earth. 
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